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Abstract 

It is shown that interference effects between velocity and density of states, which occur as elec- 
trons move along open orbits in the extended Brillouin zone, result in a change of wave functions 
dimensionality at Magic Angle (MA) directions of a magnetic field. In a particular, we demonstrate 
that these ID — » 2D dimensional crossovers result in the appearance of sharp minima in a resistiv- 
ity component p zz (H,a), perpendicular to conducting layers, which explains the main qualitative 
features of MA and Angular Magneto-Resistance Oscillations (AMRO) phenomena observed 
in low-dimensional conductors (TMTSF) 2 X, (DMET-TSeF) 2 X, and a-(BEDT-TTF) 2 MHg(SCN) 4 . 
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Low-dimensional organic conductors (TMTSF) 2 X (X = PF 6 , C10 4 , ...), (DMET-TSeF) 2 X 
(X = AuCl 2 , ...), and a-(BEDT-TTF) 2 MHg(SCN) 4 (M=K, Tl, ...) exhibit a number of 
unconventional angular magnetic oscillations [1-24] related to open quasi-one-dimensional 
(Q1D) sheets of Fermi surface (FS) in a metallic phase [1-3], 

e ± (p) = ±v F (p x TPf) - 2t b cos(p y b*) - 2t c cos(p z c*) , p F v F > t b > t c , (1) 

where +(— ) stands for the right (left) sheet of the FS; v F and p F are the Fermi velocity 
and Fermi momentum along conducting x-axis, respectively; t b and t c are the overlapping 
integrals between conducting chains; Ti = 1. Most unconventional angular oscillations in a 
metallic phase - the so-called Danner-Kang-Chaikin oscillations [17], the third angular effect 
[18-20], and the interference commensurate (IC) oscillations [20,21] - have been explained 
in term of Fermi liquid (FL) approach to anisotropic Q1D spectrum (1) (see Ref. [17], Ref. 
[25], and Refs. [26,27], correspondingly). 

On the other hand, despite the fact that all experimentally observed " magic angle" (MA) 
phenomena [5-16] and AMRO [22-24] are related to MA directions [4,28,35] of a magnetic 
field, 

tana = (n/m) (b*/c*) , H = (0, ifsina, Hcosa) , (2) 

(where n and m are integers) corresponding to periodic electron orbits in (p y ,p z )-pl&ne 
[4,35], there is no good agreement between the numerous theories of MA phenomena [28-39] 
and experiments [5-16] in a metallic phase. There exist even experimental evidences that, 
although some MA effects in a metallic phase [7,16] are of FL origin, the others [3,12-14] may 
significantly break FL picture. So far, the best qualitative agreements have been achieved 
between the prediction of Ref. [4] and the minima in onset magnetic fields for field-induced 
spin-density- wave phases observed at MA directions of the field (2) [8,37]. 

The goal of our Letter is to demonstrate that electron wave functions, corresponding 
to open FS in a realistic tight-binding model of Q1D spectrum with electron hoping only 
between the neighboring atomic sites, 

e ± (p) = ±v x (p y ) \p x TPx{Py)\ - 2t c cos(p z c*) , p x (Py) = Pf + 2t b cos(p y b*)/v F , (3) 

change their dimensionality from ID to 2D at MA directions of a magnetic field (2) with 
m — 1: 

tana = n (b*/c*) . (4) 

In particular, we show that, in the absence of Landau level quantization for open FS (3), the 
other quantum effects in a magnetic field - Bragg reflections result in ID — > 2D dimensional 
crossovers at MA directions of the field (4). 
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In other words, electron wave functions, which are localized on the conducting chains at 
arbitrary directions of a magnetic field [4,40], become 2D (i.e., localized on some planes) 
at the MA directions of a magnetic field (4). As shown below, non-trivial physical origin 
of these ID — > 2D dimensional crossovers is related to the interference effects between 
velocity component along z-axis, v z (...), perpendicular to conducting (x,y)-planes, and the 
density of states, i> x (...). These interference effects occur as electrons move along open FS 
(3) in the extended Brillouin zone and are qualitatively different from that responsible for 
IC oscillations [27,26]. Using this finding, we demonstrate that it is possible to explain 
the appearance of MA [7,13,15,16] and AMRO [22-24] minima in resistivity component 
p zz (H,a), perpendicular to conducting planes in (TMTSF) 2 X, (DMET-TSeF) 2 X, and a- 
(BEDT-TTF) 2 MHg(SCN) 4 compounds in the framework of FL approach. We also hope 
that suggested in this Letter ID — > 2D dimensional crossovers will be the key points in 
further FL and non-Fermi-liquid (n-FL) theories of more complex MA phenomena. 

At first, let us discuss how ID — > 2D dimensional crossovers can lead to the appearance of 
MA minima in p zz (H, a) using qualitative arguments. For electrons localized on conducting 
x-chains [4,40], it is natural to expect that conductivity component cr zz (H,a) is zero in 
the absence of impurities (i.e., at 1/r = 0) and decays as 1/t 2 uo 2 c (H) ~ 1/H 2 at high fields. 
[Here, uj c (H) is one of the cyclotron frequencies related to electron motion along open FS (3), 
r is an electron relaxation time.] If, at MA directions of the field (4), electron wave functions 
become delocalized, then a zz (H, a) is expected to have similarities with conductivity of free 
electrons at H — 0. Therefore, a zz (H,a) has to saturate at high magnetic fields and is 
expected to be proportional to r. Below, we demonstrate that this qualitatively different 
behavior of a zz (H,a) at MA directions (4) is indeed responsible for the appearance of MA 
minima in p zz (H,a). 

To develop a quantitative theory, we make use of the Peierls substitution method [41] for 
open electron spectrum [42,4]: p x — > —id/dx, p y — > p y — (e/c)A y , p z — > p z — (e/c)A z . 
It is convenient to chose vector potential of the magnetic field (2) in the form A = 
(0, Hxcosa, —if a; sin a), where Hamiltonian (3) in the vicinity of p x ~ pp can be ex- 
pressed as 



e+(p) = Va 
with 



Py b - 



. d 

dx 



uj b (a)x 

Pyb ~ 



V F 



-2t r cos 



p z c H 

vf 



(5) 



Ub(a) = eHvpb* cos a/c , oo c (a) = eHvpc* sin a/ c (6) 



being cyclotron frequencies of electron motion along y-axis and z-axis, respectively. 

An important difference between Hamiltonian (5) and the Hamiltonians [27,42,4] studied 
so far is that velocity component along the conducting x-chains [i.e., operator of the density 



3 



of states, v x (...)} depends on p y and x. Although in this case the operators v x {...) and d/dx 
do not commute, nevertheless one can ignore this fact if the quasi-classical parameter 

At c /uo b {a) > 1 . (7) 

It is possible to make sure that, if one represents electron wave functions in the form 

uj b {a)u 



^e{x,Py,Pz) = eX P^ P 



Py b* ~ 



V F 



du) A(x,p y ,p z ) , 



(8) 



then the solutions of the Schrodinger equation for Hamiltonian (5) can be written as 



A(x,Py,Pz) = 



\jv x [p y b* - 



exp i 



e du 



V F 



v x [p y b* 



VF J 



COS [p 2 c* 



V X [Pyb* 



LU c (a)u j 

du 



VF 



(9) 



[In Eq.(9), we normalize wave functions by the standard condition, / ip tl (x)ip t2 (x) dx = 
8(ei — €2), and make use of the inequality (7)]. 

Let us demonstrate that suggested in the Letter ID — > 2D dimensional crossovers directly 
follow from Eq.(9). It is possible to prove that in the limiting case, where v x [...} — vp = const, 
wave functions (8,9) are always localized on conducting x-chains (see Refs. [40,4]). Below, we 
show that an account of p y - and ^-dependences of the density of states, v x \p y b* — u> b (a)x / vf], 
in Eq. (9) lead to de-localization crossovers at MA directions of a magnetic field (4). For 
this purpose, we calculate z-dependence of electron wave functions at z — Nc* (where N is 
an integer) by taking a Fourier transformation of the second exponential function in Eq.(9): 



$(x,p y ,z = Nc*) = J -^exp^Nc*) exp[2it c 1 



2tt 



V X [Pyb* - 



du 



(10) 



After straightforward calculations, z-dependence of electron wave-functions (10) can be ex- 
pressed as 



$(x,p y ,z = Nc*) = exp[-i(p(x,a,p y )N} J_ N 2t c ^/i 1 2 (x, a,p y ) + I$(x,a,p y ) 
where 



(11) 



h(x = 2nM v F /uj b (a),a,p y ) = ^2 



M o r 2nv F /u> b (a) cos 



+ 2ttM 



M=0 " u Vx 

(pyb* - uj b (a) u/v F ) 



M r2-KV F /lJJ b (a) sm 

I 2 (x = 2TiM v F /uj b (a),a,p y ) = ^2 



VF ^b(ct) 



du 



du , (12) 



M=0 " u Vx 

(p y b* - u b (a) u/v F ) 

with Jn(-) being the Bessel function [43]; M is an integer. According to the Bessel 
functions theory [43], Jn(Z) is an oscillatory function of the variable N at N < \Z\, whereas 
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it decays exponentially with N at N > \Z\. Thus, one can conclude that wave functions 
(10)-(12) are extended along z-direction if at least one of the functions ij(...) in Eq.(12) is 
not restricted [i.e., if \Ii(M ,a,p y )\ — > oo as M — > oo]. In the opposite case, where both 
functions ij(...) (i = 1,2) are restricted by the conditions \Ii(M 0: a,p y )\ : \I 2 (M 0: a,p y )\ < 
Imax, electron wave functions (10)-(12) exponentially decay with the variable z at \z — 
Nc*\>2I max . 

Note that functions (12) are written in the form of summations of infinite number of 
electron waves corresponding to electron quasi-classical motion in different Brillouin zones 
in the extended Brillouin zone picture. Therefore, the physical meaning of summations 
in Eq.(12) is related to the interference effects between velocity component along z-axis, 
v z = —2t c c*sm(p z c* + uj c (a)u/vF), and the density of states, v x [p y b* — ou b (a)u/vF], which 
occur due to Bragg reflections as electron move in a magnetic field along open orbits. 
As it is seen from Eq.(12), angular dependent phase difference between electron waves, 
2irMuj c (a)/Lu b (a), is an integer number of 2ir only at MA directions (4) of a magnetic field 
(2), where uj c {a) = n cu b (a), with n being an integer. Therefore, one comes to the conclu- 
sion: at arbitrary direction of a magnetic field, the destructive interference effects in Eq.(12) 
result in exponential decay of electron wave functions (10)-(12) along z-axis, whereas, at 
MA directions, the constructive interference effects cause to derealization of wave functions 
along z-axis. 

To calculate conductivity <j zz (H,a), let us introduce the quasi-classical operator of the 
velocity component v z (...) in a magnetic field [27]: 

v z (p z ,x) = -v° z sin[p z c* + u° c (a)x/v F ] , v° z = 2t c c*. (13) 

Since wave functions (8)-(10) and the velocity operator (13) are known, one can calculate 
<t zz (H, a) by means of Kubo formalism. As a result, one obtains 



d(b*ui 



TLU b (p y + Mi, a) 



, (14) 

Py 



a zz (H,a) ~ (— / d(b*u) / \ ^exp - / 

\v x {p y ) J-oo u b (py + u,a) L Ju 

where 

u b (p y ,a) = uj b (a) [v x (p y )/v F ] , u c (p y ,a) = u c (a) [v x (p y )/v F ] , n(a) = u c (a) / 'u b (a) , (15) 

< ... > Py stands for avereging procedure over variable p y . 

After straightforward but rather complicated integrations, Eq.(14) can be rewritten as 



f U f[y + uMH)]d Ul 
Jo 



x{exp J J[y + M 1 ft fe ( J HjjdM 1 -1 

' y 

f(y) = VF/v x (y)-l, h b (H) = uJ b (a)T , h c (H) = u c (a)r . (16) 
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Since in Q1D case p zz (H, a) ~ l/a zz (H, a), Eq.(16) solves a problem to define p zz (H, a) for 
electrons with open spectrum (3) in an inclined magnetic field (2) [44]. 

To make our results more intuitive, we consider the most important limiting case - a so- 
called clean limit, where u c (a) r>l. In this case, Eq.(16) can be significantly simplified: 



(H,a) 



2 f C * A 2 n ( 

tan 2 a — > — M 
\2b*J ^ n 2 V 



1 + K«)tF V26V n 2 Vl + M«)r] 



(17) 



1 + [c<j c (q;) — nuJb(a)} 2 T 2 1 + [c<j c (q;) + nu0b(a)] 2 T 2 
where where A n are the Fourier coefficients of function f(y) = vp/v x (y) — 1: 

^v = - / +?r /(y) cos(TVy) dy . (18) 

7T J-7T 

Eq.(17) directly demonstrates MA maxima in a zz (H,a) [i.e., minima in p zz (H,a)] re- 
lated to minima of denominators which occur at u c (a) = nu b (a) [i.e., at MA directions of 
the field (4)]. In Fig.l, we present numerical simulations of Eqs.(17),(18) for three quali- 
tatively different variants of Q1D spectra (3) corresponding to (TMTSF) 2 PF 6 , a-(BEDT- 
TTF) 2 MHg(SCN) 4 , and (TMTSF) 2 C10 4 conductors. As it is seen, (TMTSF) 2 PF 6 exhibits 
only one MA minimum, whereas the last two compounds exhibit several MA minima with 
large indexes n in Eq.(4). We stress that this qualitative feature as well as a doubling of 
a period of MA minima in the case of (TMTSF) 2 C10 4 are in a good agreement with the 
existing experimental data [5,7,13,16,22]. 

We point out that the existing alternative model to describe MA and AMRO effects 
in p zz (H,a) - a so-called Osada model [30], which is very important from methodological 
and historical points of view, in our opinion, does not have a direct physical meaning. 
The reason is that the transfer integrals t„ jm in Ref.[30] are expected to be exponentially 
small in the framework of a realistic tight-binding model [1] of low-dimensional electron 
spectra. Moreover, as it follows from Eq.(17), a hypothesis [30] that it is possible to introduce 
some "effective transfer integrals", £„, m , in a linearized electron spectrum (1) and to use 
such linearized spectrum while calculating p zz (H,a) is incorrect. Indeed, weighting factors 
in Eq.(17) depend on magnetic field orientation (i.e., on tana) and, thus, their physical 
meanings are completely different from some "effective transfer integrals", t n>m , postulated 
in Ref.[30]. 

In conclusion, we hope that ID — > 2D dimensional crossovers suggested in the Letter 
will be key points in further theories describing more complicated MA phenomena. In this 
connection, we point out that there exist three main scenarios for MA phenomena: FL one 
[4,28-33] based on Gor'kov [42,45] and Chaikin [46] approach to Q1D conductors, weak n- 
FL one [28,35], where electron-electron scattering processes depend on a magnetic field, and 
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n-FL Princeton scenario [12-14,34,36], where MA direction of a magnetic field correspond 
to FL versus n-FL crossovers. Very recently, two novel exotic mechanisms [39,47] have been 
suggested to account for MA and AMRO effects. 

This work was supported in part by National Science Foundation, grant number DMR- 
0076331, the Department of Energy, grant number DE-FG02-02ER63404, and by the INTAS 
grants numbers 2001-2212 and 2001-0791. One of us (AGL) is thankful to N.N. Bagmet and 
E.V. Brusse for useful discussions. 
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FIG. 1: Resistivity component p zz (a) = l/a zz (a), perpendicular to conducting layers, calculated 
by means of Eqs.(17),(18),(5). For a-(BEDT-TTF) 2 MHg(SCN) 4 (lower curve) and (TMTSF) 2 PF 6 
(middle curve) conductors, we use model electron spectrum e(p x ,p y ,p z ) = 2t a cos(p x a/2) + 
2tb cos(p y b*) + 2t c cos(p z c*) [1] with weak, t a /tf, ~ 3, and strong, t a /tb = 8.5 [27], Q1D anisotropics, 
respectively. For (TMTSF)2C104 conductor (upper curve), we take into account anion ordering [1] 
and, thus, use the spectrum e(p x ,p y ,p z ) = 2t a cos(p x a/2) + \J[2tb cos(p y b)] 2 + A 2 + 2t c cos(p z c) [1] 
with t a /tb = 8.5 and A = 0.2tb- In all three cases, we use the value w&(0)t = 15. 
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